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Abstract. Straight-field- line coordinates are very useful for representing magnetic 
fields in toroidally confined plasmas, but fundamental problems arise regarding their 
definition in 3-D geometries because of the formation of islands and chaotic field 
regions, ie non-integrability. In Hamiltonian dynamical systems terminology these 
coordinates are similar to action-angle variables, but these are normally defined only 
for integrable systems. In order to describe 3-D magnetic field systems, a generalisation 
of this concept was proposed recently by the present authors that unified the concepts 
of ghost surfaces and quadratic-flux-minimising (QFMin) surfaces. This was based 
on a simple canonical transformation generated by a change of variable 9 = 6{Q,C), 
where 9 and £ are poloidal and toroidal angles, respectively, with 9 a new poloidal 
angle chosen to give pseudo-orbits that are a) straight when plotted in the £, plane 
and b) QFMin pseudo-orbits in the transformed coordinate. These two requirements 
ensure that the pseudo-orbits are also c) ghost pseudo-orbits. In the present paper, 
it is demonstrated that these requirements do not uniquely specify the transformation 
owing to a relabelling symmetry. A variational method of solution that removes this 
lack of uniqueness is proposed. 
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1. Introduction 

Recent calculations [lj of heat diffusion along chaotic field lines show that the isotherms 
correspond very closely with the "approximate" magnetic surfaces, associated with 
magnetic island chains, known as ghost surfaces [2J. These surfaces include the "X- 
point" and "O-point" closed field lines of their associated islands. (By "O-point" field 
line we mean either the elliptically stable field line at the center of an island or its 
hyperbolically unstable continuation if it has undergone a period-doubling bifurcation.) 
Closed field lines make the magnetic action stationary, the hyperbolic X-point field lines 
in the chaotic separatrices being minima and the O-point field lines being minimax or 
saddle points of the action. Ghost surfaces are constructed by interpolating smoothly 
between these two closed-field-line classes by evolving the O-point field lines into the 
X-point field lines along paths of steepest descent of action, thus generating a family of 
"pseudo-orbits," i.e. paths that come close to making the action stationary. 

Ghost surfaces have nice mathematical properties but are difficult to construct 
and have no obvious physical interpretation. An alternative approach to defining 
approximate magnetic surfaces passing through magnetic islands, is to use the quadratic- 
flux- minimizing (QFMin) surfaces introduced by Dewar, Hudson and Price (21 E]. These 
surfaces have the computational attraction of being easy to construct using (pseudo) 
field-line tracing methods, and the physical attraction of being defined in terms of a 
measure of the magnetic flux transport through the surface, but have been found to 
exhibit undesirable distortions in some circumstances. Thus a unified approach that 
combines the best features of ghost and QFMin surfaces is desirable. 

Both QFMin surfaces and ghost surfaces can be formulated in terms of the action 
gradient, but the action gradient is coordinate dependent. We have recently [I] exploited 
this coordinate dependence by finding the conditions under which a transformation from 
a given "old" poloidal coordinate 9 to a "new" poloidal angle O makes ghost and QFMin 
surfaces identical (a process we call reconciliation). An added benefit of this construction 
is that it makes the pseudo-orbits straight when plotted in the (, O Cartesian plane. As 
this is similar to the way the action-angle transformation for integrable systems makes 
the true orbits rectilinear we term this a generalised action-angle transformation. 

Our principal motivation is finding an optimal generalisation of straight-field-line 
magnetic coordinates in toroidal systems. As it is well known (see e.g. [5] and references 
therein) that magnetic fields can be described as 1 |-degree-of- freedom Hamiltonian 
systems we shall build our treatment upon standard classical mechanics, as for instance 
in jH], with the poloidal angle 9 as the generalised coordinate and the toroidal angle 
( as the "time". Also, in this paper we use only the Lagrangian approach to classical 
mechanics as it is somewhat simpler in the single-torus case we study. The translation 
between the Lagrangian and Hamiltonian approaches in the context of our generalised 
pseudo-orbit approach is developed in [7j. 

An analogue of our reconciliation prescription has recently been implemented [8] 
for a discrete-time dynamical system, an iterated area-preserving map (the standard 
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or Chirikov-Taylor map), as a model problem. A variational approach was used to 
perform numerical experiments aimed at finding whether the prescription can reconcile 
ghost and QFMin almost-invariant curves at arbitrary nonlinearity. ( Unreconciled ghost 
and QFMin curves have been constructed for standard-map nonlinearity parameters up 
to k = 100 [9].) Reconciliation transformations were successfully constructed for quite 
high nonlinearity (k ~ 1) after it was realised that there was a lack of uniqueness in the 
prescription of [4] that could be fixed by reducing the number of Fourier basis functions 
appropriately. 

In the present paper we identify this lack of uniqueness as due to a relabelling 
symmetry. Relabelling symmetries also occur in MHD and fluid dynamics [TO] and are 
analogous to gauge symmetries in physical field theories [TT| . Thus, making the new 
poloidal angle unique is analogous to fixing a gauge to make the representation of a field 
unique. We propose a dual-objective-functional variational method for constructing a 
unique reconciliation transformation. 

In section [2] we review the basic classical mechanics concepts required and introduce 
the concept of periodic pseudo-orbits and associated almost-invariant tori. In section [3] 
we give a concise derivation of the reconciliation conditions found in |T| and in 
section[4]we demonstrate that, given one solution satisfying the reconciliation conditions, 
an infinity of solutions may be generated by relabelling the points at which the 
pseudo-orbits cross the ( = surface of section. A primary objective functional 
that respects the relabelling symmetry and, when minimised to zero, gives valid 
reconciliation transformations is presented in section [5| Approaches for fixing the non- 
uniqueness problem are discussed in section [6] where we propose a secondary objective 
function, whose minimisation with respect to reconciliation transformations generated 
by relabelling transformations will both remove (or at least reduce) the non-uniqueness 
and ensure invertibility of the reconciliation transformation. Some areas for further 
research are briefly indicated in section [7j 

2. A simple canonical transformation 

A (p T ,q r ) -periodic path is defined in the 9,( plane by the curve 9 = subject to 
the periodicity condition + 2^q r ) = i9(C) + 2irp T (p r and q T mutually prime integers). 
Physical examples of such paths are the elliptic and hyperbolic closed field lines (periodic 
orbits) passing through the O and X points of a magnetic island chain formed through 
the resonant destruction of a rational surface with rotational transform h = p r /q r (safety 
factor q = q v /pi), but we also consider pseudo-orbits — paths that are "not quite" 
physical. 

As our theory is based on variational principles, we also consider variations 5$ 
of paths away from either physical or pseudo-orbits. For instance, the action integral 
S defined on an arbitrary (p r , g r )-periodic path is defined as a functional of the path 
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function $ by the integral 

S[&\= / L(^#(c),C)rfC, (i) 



where L = L(0, 6 1 , £) is the Lagrangian. Varying ^ in (JlJ) and integrating by parts we 
find the functional derivative of S as the coefficient of 5$ in SS, 

where Lg and L e - denote the partial derivatives of L with respect to its first and second 
arguments, respectively. In the following we refer to 5S/56 as the action gradient as 
it can be thought of as the generalisation of the gradient of a function to the infinite- 
dimensional space of path functions i?(C)- The action gradient SS/S9 can also be shown 
to play the role of a phase-space flux density (or magnetic flux density in the case of 
field-line dynamics.) 

Hamilton's Principle [6.J is the statement that S is stationary (5S/58 = 0) on 
physical orbits, i.e. the true equation of motion is obtained by setting the action gradient 
to zero. We shall term paths for which the action gradient is not zero, but in some sense 
small, pseudo- orbits. 

We also term toroidal surfaces composed of families of pseudo-orbits almost- 
invariant tori, specific cases being ghost tori when the pseudo-orbit families are 
constructed by an action-gradient flow joining true orbits, and QFMin tori when the 
pseudo-orbit families are constructed variationally to minimise the "quadratic flux" 
^JJ(5S/59) 2 d9d(, which is a measure of flux transport through the almost-invariant tori. 
The Euler-Lagrange equation for QFMin pseudo-orbits is (d/dQ8S/89 = 0, implying 
that the action gradient is constant along each QFMin pseudo-orbit (which includes the 
case of a true orbit, when the constant is zero). For details see [7j. 

The momentum canonically conjugate to 9 is 

I = L e . (3) 

Following [1] we shall seek to reconcile the QFMin and ghost formulations on a single 
almost-invariant torus by using a canonical transformation generated by a change of 
generalised coordinate from 9 to a new poloidal angle 0, 

= 0(9,0, ( 4 ) 
transforming the Lagrangian by the point transformation pp 33, 386] L new (0, 0, () = 
L \d{9,9,()- (As this is simply a statement that the Lagrangian is invariant under the 
transformation, henceforth we leave the subscripts, "old" and "new" implicit as it is 
clear from the arguments which is meant.) We assume that the mapping 9 i— y 9 is 
invertible, implying that 9(0, () is a monotonic (increasing) function of for all (. 

Equation Ml) generates a canonical transformation 9,1 i-> 0, J, where the new 
momentum conjugate to is J = dL/dQ. Differentiating (Jij) we find 

9 = (d ( + Qd e )9(Q,C)^9 (: + Q9e, (5) 
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where 9® and 9^ denote the partial derivatives of the transformation function 9 with 
respect to its first and second arguments, respectively. (Invertibility implies that 
9 e > 0.) Thus L = L(9, C + 9 e , (), giving 

J = 9 e L 9 = 9 Q I . (6) 

Note that this canonical transformation deforms the whole phase space by rescaling the 
momentum variable by an amount that varies with 9. While this is rather drastic, it 
does not matter for our current purposes as we are only interested in the neighbourhood 
of a single phase-space torus. 

As the action along any path is invariant under the transformation, the variation 
5S is also invariant, implying 5Q5S/5Q = 595S/59 for all 59. Noting that, from Q, 
59 = #e(©,£)50, we immediately see that the action gradients in the new and old 
variables are related by 

5S _ 5S 

5Q- 9e J9- (7) 
3. Reconciliation of QFMin and Ghost surfaces 

The Euler-Lagrange equation [7j deriving from the quadratic flux minimisation (QFMin) 
principle in the new coordinate [varying the family of paths in 0, ( space making up a 
trial torus, with a given transformation function 9((, 0)] is 
d 5S 

That is, the action gradient 5S/5Q is constant on each individual member of the family 
of pseudo-orbits that makes up the (p r , q T ) almost-invariant torus under consideration. 
Denoting this constant by v we have 

% = K o) , (9) 

where ©o is the initial value of on a pseudo-orbit, which we here use as a pseudo-orbit 
label. The function ^(©o) is constant along each pseudo-orbit but varies in an oscillatory 
fashion across the pseudo-orbit family, passing through zero at the action-minimizing 
and minimax true orbits [9]. 

Ghost pseudo-orbits are defined in the new coordinates by a gradient flow driven 
by the action gradient, 
DQ 5S 

DT = m' (10) 

where the evolution variable T is a label for ghost pseudo-orbits which goes from — oo 
to +oo or vice versa depending on whether the evolution is up or down the action 
gradient. (We need to consider both cases to fill in the gaps between action-minimising 
and minimax orbits, and we use the notation DT to emphasise that this flow is across 
the family of pseudo-orbits rather than along them like the pseudo-dynamics generated 
by d/d(.) 
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By reconciliation we mean that these two classes of pseudo-orbits are equivalent, 
so T and O are functionally dependent: T = T(0 O ), DT = T'(Q )DQ . Eliminating 
8S/5Q between Q and (10) we find the reconciliation condition 



^ = T'(e Me ) , (ii) 

which puts an important constraint on the reconciliation transformation it must be 
such that DQ/DQq is independent of (. (Assuming this can be satisfied, it then relates 
the ghost pseudo-orbit evolution parameter T to the QFMin pseudo-orbit label 6 .) 
But, at C = 0, DQ/DQq = 1 by definition. Thus the reconciliation condition implies 
DQ/DQo = 1 for all C 

This condition implies that 6(C|6 ) must separate in the form 6 + /(C), with / 
arbitrary except for periodicity requirements. The simplest and most natural choice to 
try is to take / linear in (, 

e(c|e ) = e +<, (12) 



with h = p/q. Equation (12) conjugates the 6 pseudo-dynamics to rigid rotation, 

0(C|0 O ) = 0(e + <, CMo = 6(Q , 0) , (13) 

in the same way that the action-angle transformation conjugates the true dynamics to 
rigid rotation [BJ, which is why we term (O, J) generalised action-angle variables, (In 
field-line terms, they might also be called "straight-pseudo- field-line coordinates".) In 
the following we propose a variational method for satisfying the new QFMin Euler- 



Lagrange equation (|8j) and the reconciliation condition in the simplified form (12), thus 
making the new QFMin torus coincide with the new ghost torus. 

4. Relabelling Transformation 



Do the reconciliation conditions rt8j) and (12) define the reconciliation transformation 



Q uniquely? We show in this section that the answer is in the negative — if there 
exists at least one solution, then there exists an infinity of different solutions generated 
by a class of transformations, 6 6, of the form 

9(6,0 = 6 + 6(6- <), (14) 

with 6(6o) any function of 6o that is 27rg r -periodic. The necessity of the 27rg r - 
periodicity restriction is to preserve 27r-periodicity in (, as can best be seen in Fourier 
representation, 

6(6, C) = 6 + ©m r sin m r (6 - <) , (15) 

mi 

(assuming odd parity, and hence a sine series) where the resonant poloidal Fourier 
indices m T are integer multiples of q r in order that corresponding resonant toroidal 
Fourier indices n r = m r h = m r p T /q r exist. 

If 6 = 6 + *C then 6 = 6 + *C> where 6 = 6(6 , 0). Thus these transformations 
rearrange the set of pseudo-orbits 6 + *C by translating them up and down in the 
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C, 6 Cartesian plane in such a way that they are all still rectilinear with slope *. As 
each Bo labels a different orbit after the transformation than it did before, we call 
transformations of the form ( 14 ) relabelling transformations. 

Composing relabelling transformations with the reconciliation transformation ^ 
gives what we now proceed to show to be an equivalence class of reconciliation 
transformations 8(0, C), 

0(6(8,0,0 = mo- (is) 



The relabelled conjugacy equation, analogous to (13), is 

0(C|0„) = 0(Q + <, C)Jo = 6(Qo, 0) • (17) 

the velocity d' and acceleration i?" transform similarly. Thus, the "unreconciled" action 
gradient §2§ is invariant under relabelling, 



6S_ 



(^>C) = ^W,0- 



Differentiating both sides of (16) with respect to gives 
6 e (6,O0e(6,O 



(19) 



0(6,0, 

6(6 - <). Using (B and (fl9[ we find the relabelling 



(20) 



where, from (16), 6e(6,0 = 1 

transformation condition for the "reconciled" action gradient (J7|), 

= I- 

Taking the total derivative of both sides of pOj ) with respect to £ and observing 
that c?9e(6o)/^C — commutes with d/d(, if @ is satisfied before the relabelling 
transformation, it will be satisfied after, and thus the relabelling transformations (14) 
generate equivalent reconciled solutions. 

Finally, consider the case of infinitesimal relabelling transformations, 

6(0,0 = 6 + 56(0- <). (21) 



Expanding (16) to first order we find the general relabelling variation 

S9(Q, = -58(6- <) Me, 0- 



(22) 



5. Variational Formulation 



We can impose condition (12) simply by using it as an ansatz for the pseudo-orbit paths 
in 6,C space [which then defines the paths in 9, £ space via Q: $(0 = #(60 + *CjC)]> 
but imposing (|8| is more difficult as it is nonlinear. In j3] we used perturbation theory, 
but this is limited to small departures from integrability. Instead we now introduce a 
variational approach that can be used to find numerical solutions for arbitrarily nonlinear 
problems by building on standard optimisation methods. 
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As the (primary) objective functional we take 



F[0\ = 



27r /qi 



q v dO( 
2vr 



2irq r 



2-nq T 



1 



d~cse 



(23) 



J tf=0(e,c),e=e o +*< 



2tt 



e 



0fi 



55 



-i 2 



where we have obtained the second form of F by changing variables from (, © to 
C,©, so that, in the first form, d/d( denotes with Go fixed, while in the second 
form d/d( = + td®. The details of the reshuffling of the limits of the integrals are 
similar to those spelt out in eq. (45)]. Also, in the second form we have used (uh 
to express SS/SO in terms of 6q(0, £) and 5S/56, which is given explicitly in terms of 
the Lagrangian L in (J5]) . The weight factor 1/Qq has been inserted to make F precisely 
invariant with respect to relabelling transformations whether or not ^ is satisfied, as 
can be demonstrating by making a change of variable from O to and using the results 
of section HI 

The objective functional F is to be minimized over all functions 6(0, () such 
that 8(0, Q is a 27r-periodic function of (, is monotonically increasing in O, and 
8(0 + 2tt,() = 6(Q,() + 27T. (Or, equivalently, such that 6q is 27r-periodic in both 
and £, its ©-average is unity, and its ^-average is zero.) Clearly F[6] > 0, with 
equality applying iff pj) is satisfied. 



6. Constraining the relabelling symmetry 

If F = for some 6(0, an infinity of new solutions can be generated by applying 



finite relabelling transformations (16) — the solution to the QFMin condition Q) is not 
unique. Without a further constraint to fix the solution, no optimisation algorithm for 
minimising F can ever converge. 

One such constraint method derives from the Fourier form of the relabelling 



transformation, (15), where it is seen that the relabelling symmetry gives us precisely 
enough freedom to constrain the resonant Fourier coefficients in the expansion of 6 to any 
desired value, thus removing the non-uniqueness in the reconciliation transformation. 

This constraint approach was implicitly adopted in the perturbation method given 
in [7], where we set the resonant Fourier coefficients in the expansion of 6 to zero at 
first and second order in nonlinearity. An analogous zero-resonant-coefficients constraint 
method was used in a nonlinear numerical study [S] , performed using the analogue of our 
primary objective functional F for the standard map. While this study gave convincing 
evidence that reconciliation transformations exist for systems with quite large islands 
and chaotic regions (standard map nonlinearity parameter k ~ 1) it was found that the 
method broke down for stronger nonlinearity because the transformation O y 6 became 
non-invertible for k 3> 1. 

As there is no compelling reason to set the resonant coefficients to zero, we propose 
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that a better approach would be to determine them by minimising a secondary objective 
function G, not invariant under relabelling, over the equivalence class generated by 
the relabelling transformations (14). As the failure of the zero-resonant-coefficients 
constraint method manifested itself by 6@ going negative, a natural choice for secondary 
objective functional is the average of the weight function 1/9q, 



which maintains monotonicity because G diverges toward +oo if 8q — > +0 anywhere in 
the integration domain. Using the Schwarz inequality for the inner product between the 
two functions a/#6 and 1 / \/~9e we can a l so show that G is bounded below by unity, so 
1< G < oo. 



Preliminary numerical experiments on minimising the analogue of (24) for the 
standard map to fix the resonant 8 m indicate that breakdown of invertibility can 
be avoided for higher k, but more systematic studies need to be performed, using 
appropriate dual-objective numerical optimisation methods, before one can conclude 
that reconciliation can be performed for arbitrary nonlinearity. 

Unfortunately, despite the ubiquity of Lie symmetries in physics, there appear 
to be no algorithms in the standard numerical optimisation texts appropriate to this 
problem. Multi-objective optimisation problems are well known, but these involve a 
trade-off between competing objectives (the Pareto problem), whereas we wish to give 
100% Pareto weight to F and minimise G only over the subspace of directions where 
F does not change. [Note that this nullspace is not precisely the same as the subspace 
spanned by the resonant Fourier modes because of the factor 6® in (22).] 

A similar relabelling problem arises in three-dimensional numerical MHD 
equilibrium calculations, where the objective function is the total plasma and magnetic 
energy and the relabelling symmetry arises from the arbitrariness of choosing the 
poloidal angle within magnetic surfaces. A numerical method [121 H3] for fixing the 
poloidal angle has been implemented in the VMEC code [Ej, using a measure of the 
width of the Fourier spectrum as the secondary objective function. Adaptation of 
this method, and other optimisation methods, to the present problem will be reported 
elsewhere. 



7. Conclusion 

We have reviewed the motivation and formulation of a recently published [5] 
unification (reconciliation) of ghost and quadratic-flux-minimizing (QFMin) surfaces by 
transforming to a new poloidal angle, and have identified a relabelling symmetry that 
makes the reconciliation transformation non-unique. We have proposed a variational 
approach using a primary objective function to satisfy the reconciliation conditions 
and a secondary objective function to fix the relabelling symmetry, giving an explicit 
expression for the gradient of the primary objective function and identifying the 
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nullspace of its Hessian operator (the space spanned by infinitesimal relabelling 
transformations) . 

Numerical validation of the method at high nonlinearity remains for further 
work. Also, the localised variational approach presented here is based on a canonical 
transformation tailored to reconciling the ghost and QFMin approaches on a single 
resonant surface using a Lagrangian approach (a point transformation). To find a global 
pseudo-magnetic coordinate system [TS1 HB] , we need to generalise this approach to find 
a canonical transformation that allows a simultaneous multi-surface optimisation. We 
anticipate that the generality of Hamiltonian methods for constructing ghost and QFMin 
surfaces [7] is better adapted to this purpose than the simple Lagrangian approach used 
in the present paper. 
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